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Abstract We discuss constraints on the effective number of neutrino species N v from recent cosmolog- 
ical observations such as CMB, LSS, BBN, including our own analysis which uses the WMAP and the 
Luminous Red Galaxy power spectrum data. We also discuss their implications on some non-standard 
cosmological scenarios such as the low (MeV-scale) reheating temperature scenario and the scenario with 
decaying particles between BBN and structure formation. 



1 Introduction 

We know very well about the energy density of relic photons from the hot Big Bang 
since the COBE measurement confirmed a Planck distribution of the background photons 
at the temperature T = 2.725 K [HE]- In the standard cosmology, for relic radiation 
(relativistic) component other than photons, we assume three generations of neutrinos 
with a Fermi distribution, but we have not directly measured them. Also, there may be 
other contribution which cannot be directly detected. Such uncertainty in our knowledge 
of the radiation energy density in the universe is often expressed as the effective number 
of neutrino species N v . N v is the radiation energy density (other than that of photons) 
normalized by the neutrino energy density for one species in the standard cosmologyo 
We can think of many candidates in particle physics/cosmology that could modify the 
standard value of N v = 3. There may be sterile neutrinos, gravitational waves, axions, 
majorons, large lepton asymmetries, MeV-scale reheating and so on. Therefore it is of 
great importance to constrain N v from the cosmological observations to probe a certain 
class of models in the particle physics/cosmology. 

Recent precise observations of the cosmic microwave background (CMB) anisotropics 
and large scale structure (LSS) make it possible to measure N u through its effects on the 
growth of cosmological perturbations. These effects come from the fact that the den- 
sity perturbation does not grow (the gravitational potential decays) during the radiation- 
dominated era. Specifically, a more relativistic degree of freedom causes greater early 
integrated Sachs- Wolfe effect on the CMB power spectrum, which leads to higher first 



1 More precisely, — (p rc \ — p-() / Pv.thm, where p 7 is the photon energy density, p le \ is the total energy 
density of photons, three active species of neutrinos and extra relativistic contribution, and p^.thm is defined 
as p^thm = (77r 2 /120)(4/ll) 4 / 3 T 4 using the photon temperature T after the electron-positron annihilation. 
PiAthm corresponds to the energy density of a single species of neutrino, assuming that neutrinos are 
completely decoupled from the electromagnetic plasma before the electron-positron annihilation takes 
place and they obey a Fermi-Dirac distribution. 
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Table 1: Comparison of N v constraints using various data set combinations. "All" refers 
to WMAP3 + other CMB + Lya + galaxy power spectrum (SDSS main sample + 2dF) 
+ SDSS baryon acoustic oscillation (BAO) + Supernovae la (SN). See Ref. [I] for details. 
SDSS (main) and Lya favor N u > 3. 

peak height. Also, since it delays the epoch of the matter-radiation equality and makes 
the horizon at that time larger, the turnover position of the matter power spectrum is 
shifted to larger scales and the power at smaller scales is suppressed. Therefore, by ob- 
serving the CMB and LSS, we can measure during structure formation. In detail, 
assuming the smallest scale relevant to our observations to be about 5Mpc, since struc- 
ture formation of that scale begins around the temperature T « 20 eV (at which the scale 
enters the horizon), these observations probe N v at T < 20 eV (if 1 Mpc is taken to be 
the smallest observable scale, as is the case with Lyman-a forest measurements (Lya) 
, it would be T < 100 eV). Thus, the CMB and LSS can measure N u independently of 
another well-known probe, big bang nucleosynthesis (BBN), which measures it in much 
earlier universe around T = O(MeV). 

There are now many works to constrain N v from cosmological observations and, sur- 
prisingly, some data sets seem to favor N u > 3 as found in Refs. [3] EJ [5] [6] (we show 
the result of Ref. pE] in Tabled]). The galaxy power spectrum of 2dF |7J favors standard 
N u = 3 but the SDSS main galaxy power spectrum [8] and the Lya data prefer N v > 3. 
Thus, we would like to cross-check these disagreements by using the SDSS luminous red 
galaxy (LRG) power spectrum [9] which has more statistical constraining power than 2dF 
or SDSS main sample (the effective volume of the LRG survey is about 6 times larger than 
that of the SDSS main galaxy sample and over 10 times larger than that of the 2dFGRS 
[9j)@ 

We describe our analysis and result in the next section. The implications for some 
non-standard cosmological scenarios are discussed in Sec. [3J We summarize in Sec. HI 

2 Recently, Ref. [10] argues that the discrepancies are due to systematic effects of scale-dependent 
biasing in the galaxy power spectrum. 
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Figure 1: A% 2 as functions of N v . The red solid line shows the result of using the WMAP 
three-year data alone and the green dashed line that of using the WMAP three-year data 
and the SDSS LRG power spectrum [TT] . 



2 Constraint on N v from WMAP and the SDSS LRG 
power spectrum 

We summarize the analysis and result obtained in Ref. [TT]. We constrain N v in the 
flat ACDM universe with the initial perturbation power spectrum which is adiabatic and 
described by a power law. This model has six cosmological parameters: the baryon density 
ujf,, the matter density uj m , the normalized Hubble constant h, the reionization optical 
depth r, the scalar spectral index of primordial perturbation power spectrum n s and its 
amplitude A (u = Qh 2 , where Q is the energy density normalized by the critical density). 
Theoretical CMB and matter power spectra are calculated by the CMBFAST code [12] 
and x 2 by the likelihood codes of the WMAP three-year data [131 1131 H5j and of the SDSS 
LRG power spectrum data [9]. We apply modelling of nonlinearity and scale-dependent 
bias as in Ref. [9] to the linear matter power spectrum before fitting to the LRG data. This 
modeling has two parameters, the galaxy bias factor b and nonlinear correction factor Q n \. 
Specifically, we connect the linear matter power spectrum P\\ n {k) and the galaxy power 
spectrum P m \{k) by P ga i(k) = b 2 {(l + Q nX k 2 )/{\ + lAk)}P lin (k). We calculate the \ 2 as 
functions of N v by marginalizing over the above parameters (6 parameters for WMAP alone 
and 8 for WMAP+SDSS). The marginalization is carried out by the Brent minimization 
[T6] modified to be applicable to multi-dimension parameter space as described in Ref. [TT] . 

We show the results of x 2 minimization in Fig. [TJ We have checked that the results 
for the standard three neutrino species agree with the WMAP [3J and SDSS [H] groups' 
analyses. For the WMAP three-year alone case, it has been checked in Ref. [TB] that the 
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best fit x 2 an d parameters agree. With regard to the WMAP and LRG combined analysis, 
our best fit parameter values for the three neutrino species are w& = 0.0222 ± 0.0007, 
oo m = 0.1288 ± 0.0044, h = 0.718 ± 0.018, r = 0.088 ± 0.029, n s = 0.958 ± 0.016, a s = 
0.770 ± 0.033 (we report here a& instead of A to compare with Ref. [9]), b = 1.877 ± 0.065 
and Q n \ = 30.4 ± 3.5. The central values fall well within the la ranges of the constraints 
derived in Ref. [9] and the la errors are almost identical to those quoted in Ref. [9]. 

The limits corresponding to A% 2 = 4 are N u < 25 for WMAP three-year data alone 
and 0.8 < N v < 8.0 for WMAP and SDSS LRG combined. Since the \ 2 functions show 
some asymmetric features, we derive 95% confidence limits by integrating the likelihood 
functions C = exp(-A X 2 /2). This yields a 95% C.L. bound of N u < 42 for WMAP alone 
and 0.9 < N v < 8.2 for WMAP+LRG. 

Our constraint on N v from the WMAP three-year CMB power spectrum and the SDSS 
LRG power spectrum, can be summarized as 



at the 95% C.L. The minimum x 2 is given at the value quite close to the standard model 
value but the upper bound is not so stringent that it is compatible with the Lya's high 
value. Since it is very difficult to remove systematic errors from the galaxy clustering and 
Lya analyses, we may have to wait for next generation CMB experiments to significantly 
improve our knowledge on N v (the PLANCK sensitivity for N u is forecasted to be ~ 0.2; 
see e.g. Ref. p2]). 

3 Implications for non-standard cosmology 
3.1 Low reheating temperature 

In this scenario, the effective number of neutrinos N v can be less than the standard value 
since neutrinos are not thermalized. We briefly explain this scenario and how N v and the 
reheating temperature Tr are related. For more details, we refer to Ref. [20] . 




The standard big bang model assumes that the universe was once dominated by thermal 
radiation composed of photons, electrons, neutrinos, and their antiparticles. The reheating 
temperature is the temperature at which the universe becomes such a radiation-dominated 
state and it is usually assumed to be so high that every particle species is in thermal 
equilibrium. In particular, neutrinos are considered to obey a Fermi distribution. But if the 
reheating temperature is as low as a few MeV, their distribution function deviates from the 
thermal one. In contrast to electrons that are always (at least until the temperature drops 
below a few eV) in thermal contact with photons via electromagnetic forces, neutrinos 
interact with electrons and themselves only through the weak interaction. The decoupling 
temperature of the neutrinos should be around 3 MeV for the electron neutrinos and 5 
MeV for the muon and tau neutrinos, respectively (the difference comes from the fact 
that the electron neutrinos have additional charged current interaction with electrons). 



0.9 < iV^ < 8.2 or N u = 3.1 



(1) 
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T R (MeV) 

Figure 2: The relation between the effective neutrino number N v and the reheating tem- 
perature T R [2D] . 

Therefore the neutrinos might not be fully thermalized and lead to N„ < 3 if the reheating 
temperature is in the MeV range. In fact, such a low reheating temperature can be found 
in many cosmological scenarios especially in order to dilute the unwanted relics such as 
the gravitinos. 

In Ref. [20], we have calculated how much neutrinos are thermalized when Tr = 
O(MeV) and have derived the relation between Tr and N v which is shown in Fig. [2]|2| 
We can convert our constraint on N u , Eq. ((T]), into the lower bound on T R using Fig. [2j 

T R > 2 MeV. (2) 



3.2 Decaying particles between BBN and structure formation 

As mentioned at the end of Sec. [2J our analysis using the SDSS LRG power spectrum 
cannot exclude N v = 4 ~ 5 which is favored by the Lya data. On the other hand, the 
BBN constraint from the recent measurements of the 4 He abundance, Y p , seems to favor 
N v ~ 3 as we summarize in Table [2j If these Lya and Y p analyses are correct, can we 
obtain a consistent cosmological scenario ? 

Recall that they measure N u at different cosmological epochs as noted in Sec. [U Thus, 
we can do this by increasing N v after BBN but before structure formation begins. This 

3 As discovered in Ref. [5D], it should be noted that the 4 He abundance, Y p , increases while N v decreases 
in this scenario. This is in contrast to the conventional non-standard N v scenario where decreasing 7V„ 
accompanies decreasing Y p . The difference occurs because non-thermal distribution of (electron- type) 
neutrinos modifies the neutron-proton conversion rate, which significantly affects Y p . 
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Table 2: Comparison of N u constraints from recent Y p measurements. We also used the 
observed deuterium abundance D/H = (2.82 ± 0.27) x 10~ 5 [25] and the BBN fitting 
formula in Ref. [26] . N v > 4 is not favored by the three recent measurements. 

corresponds to the temperature between about 1 MeV and 100 eV or the time between 
Is and 10 8 s. In Ref. [27], we have investigated models with decaying particles which 
have a lifetime of O(l-10 8 ) s and can contribute to the required additional N u . We have 
considered (a) saxion decay into two axions; (b) gravitino decay into axino and axion; (c) 
Dirac right-handed sneutrino decay into gravitino and right-handed neutrino, and have 
shown that there are parameter regions which can make AN U ~ 1 and do not affect other 
observations. 

Finally, we would like to comment that the reason why the Lya data favor N u > 3 is 
connected to the <r 8 value measured by the Lya which is higher than the WMAP3 value, 
as is pointed out by Ref. [4J. In the standard ACDM model (i.e. N u = 3), WMAP3 
alone gives <7g = 0.76 ± 0.05 [3] while combining it with the Lya gives much higher value 
cx 8 = 0.85 ± 0.02 [4]. However, if we remove the assumption N u = 3 and increase N v , the 
best fit value of <r 8 from WMAP increases (we can see this explicitly in Fig. [3]) and more 
agrees with the Lya data. 

Now, turning the argument around, we can say that observations which give relatively 
higher er 8 basically favor N u > 3. We know that some of the weak lensing observations 
indicate high cx 8 . For example, Ref. [28] gives <r 8 = 0.84 ± 0.07 using a ground-based 
telescope and Ref. [29] gives a 8 = 0.95^^1 using a space telescope. Ref. [30] combined the 
latter with the Lya data and obtained, naturally, a high value a 8 = 0.87 ±0.05. Although 
the statistical significance is not very high, the weak lensing data too give motivation for 
the scenario presented in this section. 

4 Summary 

We have discussed the cosmological constraint on the effective neutrino species N u and its 
implications for some of non-standard cosmological scenarios. Our analysis showed that 
the constraint from the CMB data of WMAP3 combined with the galaxy power spectrum 
of the SDSS LRG sample is 0.9 < N u < 8.2 at 95% C.L. The lower bound on N v can be 
converted to the lower bound on the reheating temperature, Tr > 2 MeV. The upper bound 
turned out to be not very stringent so the preference of N v > 3 indicated in the Lya data 
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Figure 3: The values of Ug when we minimize \ 2 f° r the WMAP three-year data alone (the 
red solid line) or together with the SDSS LRG power spectrum (the green dashed line). 

cannot be ruled out. This motivates us to investigate a scenario in which N u = 3 during 
BBN (as constrained from the 4 He abundance) and N v > 3 before structure formation. 
N v > 3 is considered to be also favored by some of the weak lensing observations which give 
high erg values as the Lya data. We found that we can increase N v by decaying particles 
within the framework of supersymmetric extensions of standard model of particle physics. 
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